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We use small angle neutron scattering to study the superconducting vortex lattice in 
La2-^Sr :E Cu04 as a function of doping and magnetic field. We show that near optimally dop- 
ing the vortex lattice coordination and the superconducting coherence length £ are controlled by a 
van-Hove singularity crossing the Fermi level near the Brillouin zone boundary. The vortex lattice 
properties change dramatically as a spin-density-wave instability is approached upon underdoping. 
The Bragg glass paradigm provides a good description of this regime and suggests that SDW order 
acts as a novel source of disorder on the vortex lattice. 



I. INTRODUCTION 

A commonality across the borocarbides,— ~— 
cupratesj^ ferro-pnictides^ heavy- fermion 7 - and or- 
ganic superconductors^ is the coexistence of magnetism 
and superconductivity. The corresponding order param- 
eters typically compete^ and often a small perturbation 
is sufficient to tip the balance between the two. For 
example, the magnetism carried by the rare earth ions 
R in the borocarbides RM2B2C can lead to nearly 
reentrant superconducting phase diagrams, and sponta- 
neously forming superconducting vortices at zero applied 
fields— Equally, vortices induced under applied fields 
may permit enhanced magnetic correlations in the core 
regions where the superconducting order parameter is 
suppressed.— ~— This idea was put forward to explain 
field-induced and -enhanced magnetic correlations 
observed in the cuprate superconductor La2- x Sr^Cu04 
(LSCO).- 5 -^ 

Although the effect of static magnetism on moving 
vortices was recently considered theoretically— little is 
known about how the presence of magnetic correlations 
affects the arrangement of vortices. Here we address 
the problem from an experimental point of view. When 
magnetism and superconductivity coexist there are at 
least four relevant length scales: the penetration depth 
A, the vortex core size £, the vortex spacing ao, and 
the magnetic correlation length £. The vortex spacing 
ao oc H~ 5 scales with the applied magnetic field n^H 
and in LSCO the magnetic correlation length can be 
tuned by varying the doping concentration. Using small 
angle neutron scattering (SANS) we have studied two dif- 
ferent regimes (see Fig. 1); (i) far away from the magnetic 
ordering where £, £ <C ag and (ii) entering the magnetic 



phase where £ ~ ao. In the first regime, where static 
magnetism is absent, the vortex lattice (VL) structure 
and core size are understood from pure fcrmiological con- 
siderations. In the second regime with static long-range 
magnetism, the vortex arrangement exhibits increasing 
disorder. We find these regimes to be well-described 
within the topical 'Bragg glass' paradigm, where disor- 
der results in an algebraic decay of the translational or- 
der of the vortices . 16 ' 17 VL disordering is usually driven 
by effects extrinsic to superconductivity such as rare 
earth magnetism in RNi2B2Cji — ii2r— or sample impuri- 
ties and crystalline defects^ In contrast, magnetic and 
SC order parameters are intertwined in LSCO; we show 
that this provides a novel and tunable source of VL dis- 
order. 

As shown in Fig. 1, the appearance of magnetism is, in 
essence, concomitant with the suppression of SANS in- 
tensity with field and underdoping. Drawing upon results 
from the literatur e) 23 ' 24 ' 29 and new observations reported 
herein, we are also able to plot the VL structure at low 
temperature T versus magnetic field H and doping x. 

II. EXPERIMENTAL METHODS 

Single crystals of La2_ a; Sr a ;Cu04 with x = 0.105- 
0.22 were grown by the traveling solvent floating zone 
method^ The static and dynamic magnetic properties 
of the samples were characterized using both neutron 
diffraction and neutron spectroscopy and good agreement 
was found between our dat a 25 ' 31 ' 32 and previously pub- 
lished results ^i^ 6 -^— 

The SANS experiments reported here were carried out 
over a series of experiments using the SANS-I instrument 
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FIG. 1. (Color online) Phase diagram of (a) the (T < 
3 K) vortex lattice structure and (b) the magnetism in 
La2- a; Sr :E Cu04, both revealed by neutron diffraction. Data 
points in (a) for x = 0.12, 0.145, 0.16 and 0.22 are from this 
work whereas x = 0.105, 0.17, and 0.20 are from Ref.lH^. 
Circular points are defined by the onset of a square VL coor- 
dination. The field-doping plane in (b), adapted from Ref . |25|. 
shows schematically the ordered SDW-moment normalized to 
that at the 1 / 8-doping. Field-induced order was first rep orted 
by Khaykovich et al^ and later confirmed in Ref. 12a . The 
temperature-doping plane shows the superconducting dome 
together with the onset of static incommensurate SDW order 
Tsdw as seen by neutron diffraction.— 1 ^ We remark that a 
similar phase diagram was proposed for YBa2Cu30 a .— The 
dashed lines indicate the samples studied in this paper. 



at SINQ^ 6 . and the Dll and D22 instruments at ILL. In 
all experiments we adopted the experimental geometry 
where external magnetic fields [IqH are applied paral- 
lel to the crystal c-axis, and almost parallel to the neu- 
tron beam. The scattered neutrons are recorded using a 
position-sensitive-detector placed behind the sample. For 
each doping, up to three different neutron wavelengths 
spanning the range X n = 5-16 A were used in order to 
cover the applied field figH range of 0.03-10 T. In all 
cases, a zero-field cooled background was subtracted from 
the field-cooled data in order to leave just the VL signal. 

Our experimental setup is shown in Fig. [2] where we 
also illustrate the relationship between the VL in the 
sample and the quantities extracted at the position- 
sensitive-detector. Typically, to observe the signal due to 




FIG. 2. A schematic diagram illustrating the experimen- 
tal geometry chosen for our experiments. A square VL in real 
space forms a two-dimensional reciprocal VL, the properties of 
which are recorded by SANS using a position-sensitive detec- 
tor. The Bragg spots in reciprocal space exhibit finite widths 
we, w q and w± that are dependent on both the instrumen- 
tal resolution and properties of the VL. These three lengths 
are estimated at the detector by recording the angular widths 
r r and ta within the detector plane, and t^ perpendicular 
to the detector plane. t u is determined experimentally by 
recording the rocking curve, and corresponds to the rocking 
curve width. 



the VL, the sample and cryomagnet are rotated together 
by angles (such as that shown by uj in Fig. 2) in order 
to bring a reciprocal lattice vector onto the Bragg con- 
dition at the detector. Due to both the finite resolution 
of the instrument, and the mosaic spread or imperfec- 
tion of the VL, the Bragg spots occupy a finite volume 
in reciprocal space, and can be described, in a first ap- 
proximation, by three widths wi, w q and w± as shown 
in Fig. 2. These three lengths correspond to the angu- 
lar widths r w , r r and ta respectively measured in the 
SANS experiment. r r and ta describe the finite size of 
the Bragg spot on the detector plane, and t w is the width 
of a rocking curve, measured by recording the Bragg spot 
intensity as a function of rotation angle uj. In the geom- 
etry shown in Fig. 2, ta (and w±) provide a measure of 
the VL orientational order about the field axis, while T r 
(and w q ) is dominated by the SANS instrumental reso- 
lution function. The contribution due to the resolution 
function is smallest for t u (and wg), which is sensitive to 
the VL correlation length along the field direction (vortex 
'straightness'). In each of our measuremenets, the sample 
and the field were rotated together through a wide range 
of rocking angles uj spanning ± 5° about the neutron 
beam. Note that this angular rotation range is typically 
narrower than the width of rocking curve in our sam- 
ples. However, by summing over the measured rocking 
angles, diffraction patterns such as those shown in Fig. 
3(a) and (b) are obtained, allowing a determination of 
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FIG. 3. (Color online) (a)-(b) Vortex lattice diffraction pat- 
terns, recorded using a position sensitive detector (see Fig. 
2), of LSCO x = 0.145 under applied magnetic fields of 0.2 T 
and 1.3 T, respectively. Notice that the Cu-O axes are along 
the diagonal, (c)-(d) Azimuthally-averaged momentum In- 
dependence of the scattered intensity, summed over rocking 
angles, for dopings x = 0.145 and 0.12, and applied magnetic 
fields as indicated. For visibility, the 0.2 T data in (c) have 
been divided by a factor of ten. Inset of (c) is a zoom on 
the 0.2 T data. The vertical bars above the diffraction peaks 
indicate the expected positions for square and regular hexag- 
onal VL coordinations. The solid black lines in (c)-(d) are 
Lorentzian fits to the data. 



the VL properties. 



III. RESULTS AND DISCUSSION 



A. Vortex lattice morphology 



Our observations of the VL structure and coordination 
can be quantified in terms of a dimensionlcss parameter 
a = 47r 2 /ioi?/($o|G| 2 ) where |G| is the magnitude of 
the reciprocal VL vector. For a regular hexagonal VL 
coordination a = \/3/2, while for a square coordination 
(7 = 1. The definition of a is useful because it does not 
require details of the positions G of Bragg peaks; only 
the magnitude |G| is needed. At low fields, the VL is 
susceptible to oricntational disorder due to impurities or 
defects in the sample, but |G| can still be measured. For 
example, in Fig. [3ja) the diffracted intensity measured 
in LSCO x = 0.145 at 0.2 T does not show well-defined 
Bragg spots. This indicates a large ta and poor orienta- 
tional order of the VL about the c-axis. Nevertheless we 
can determine the VL coordination by averaging over the 
azimuthal angle ip so that the diffracted intensity /(|q|) 
becomes a function of |q| only, see Fig. (3Jc). |G| is de- 
termined from the peak position in the |q|-dcpcndence. 
Fitting a Lorentzian lineshape to these 0.2 T data, yields 
cr = 0.88(2) very close to the value expected for a hexag- 
onal VL. At hqH = 1.3 T in Fig. EJTj), we find four 
Bragg spots with G along the Cu-0 bond directions and 
(7 = 0.99(1), indicating not only an improved VL orien- 
tational order but moreover a square coordination. As 
shown quantitatively in Fig. 4(a), the VL coordination 
in underdoped LSCO x = 0.145 changes steadily from 
hexagonal to square over the range hqH = 0.2 to 0.8 T. 
In contrast, on the optimally- and over-doped side of the 
phase diagram (x > 0.17), the VL structure becomes 
square by fi H « 0.4 T (Fig. Ufa)). 

A minimum in either the Fermi velocity vp(k) or 
the superconducting gap A(k) are well-known sources of 
field-driven hexagonal-to-square VL transitions^ - — In 
both LSCO and YBa 2 Cu 3 O y (YBCO) the band structure 
is predominantly two-dimensional but with some c-axis 
dispersion near the k = (it, 0)-point . 43 ' 44 It was previ- 
ously suggested that the VL morphology should be un- 
derstood from the Fermi surface topology near the (7r,0)- 
pointi2S A crucial difference between the band struc- 
tures of LSCO and YBCO is that in LSCO, a van Hove 
singularity crosses the Fermi level (cf) at (it,0) some- 
where between x = 0.17 and 0.22 t 44 i 45 leading to a 
huge t>i?-anisotropy. This offers an explanation as to 
why the square VL is oriented along the (7T,0) direc- 
tion in LSCO as opposed to the nodal (7r,7r)-direction in 
YBCO^-ii^-iS As LSCO is underdoped, the van Hove 
singularity is pushed further away from e^,— leading to 
a smaller vp -anisotropy and consequently a larger field is 
required to form the square VL, as indeed observed for 
LSCO x = 0.145 (Fig. 4a). 
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FIG. 5. (Color online) Superconducting coherence length 
£ in LSCO extracted from magneto-resistance (solid blue 
points)^ specific heat (open black points)^ and our SANS 
measurements (solid red points) and plotted as a function of 
the hole concentration x. The high-field magneto-resistance 
study measures the upper critical field H C 2 at low tempera- 
tures and we used H C 2 = 3>o/ (27r£ 2 ) to estimate the coherence 
length £. On the other hand, the specific heat and the SANS 
experiments were carried out at fields smaller or compara- 
ble to H C 2- To extract the superconducting coherence length 
from the SANS data we used the Clem model for the VL form 
factor. Notice that in presence of vortex lattice disorder, the 
Clem model will overestimate the coherence length (see text). 
This may explain why the coherence length extracted from 
the SANS data lies systematically above the specific heat and 
magneto-resistance measurements. All lines are guides to the 
eye. 



FIG. 4. (Color online) (a) Dimensionless constant a, de- 
fined in the text, as a function of magnetic field for LSCO 
x = 0.145, 0.16 0.17 (Ref. [H), 0.20 (Ref. [H), and 0.23. 
a = \/3/2 is expected for a hexagonal lattice and a square 
vortex lattice has a = 1. The change form a — y/3/2 to 1 
therefore reveals the a hexagonal-to-square transition of the 
vortex lattice structure (see also Fig. 1). (b) SANS intensity 
I versus applied magnetic field for NCCO x — 0.15 (Ref. l37f ) 
and La 2 - ;c Sr ;c Cu04 with x = 0.105 (Ref. S3), and 0.145-0.23, 
(this work). For clarity, the intensities for each of the compo- 
sitions have been given an arbitrary vertical offset. Solid lines 
are fits to the Clem model form factor where the supercon- 
ducting coherence length £ is the only parameter (see Fig. 5). 
Dashed lines indicate power law dependencies; / ~ H~ 5 and 
I ~ H~ 2 . Notice that the Bragg glass paradigm for vortices 
in the presence of disorder is consistent with a cross over from 
I ~ H-° b to I ~ H~ 2 , see text. 



B. Diffracted SANS intensity 

We define the VL intensity of the first order diffraction 
peak / as the sum of the area under 7(|q|), which is itself 
a sum over rocking angles. Overlap measurements of / vs 



HoH, shown in Fig. 4(b) for x = 0.105-0.23, were done 
whenever the neutron wavelength A„ was changed. For 
x > 0.16, intensity could be observed up to the highest 
applied field 10 T. By contrast, for x = 0.145 no intensity 
was observed above the quantum critical field (poH = 7± 
1 T in our sample) for SDW order i 25 ' 26 i 31 It is, however, 
still possible that the VL extends slightly into the SDW 
ordered phase. This is the case in LSCO x — 0.105, 22 
where a 3D vortex lattice exhibits / oc H~ 2 (Fig. Hfb)) 
over two decades of intensity and co-exists with short 
range SDW order— at very low fields H c \ < HqH < 
0.2 T. VL intensity is also observed in LSCO x = 0.12; 
a compound where long-range SDW order exists already 
in zero fields At fioH = 0.05 T, the |q|-depcndence 
of the intensity /(|q|) (Fig. Eld)) suggests a hexagonal 
VL coordination. To the best of our knowledge, this 
provides the first evidence by SANS of a VL co-existing 
with SDW-order at the 1/8-anomaly. The field range 
of co-existence is small; on increasing the field to just 
(j,qH = 0.1 T, Fig. [3jd) shows that the VL signal has 
already fallen to the background level. Notice that the 
fields 0.1-0.2 T are much smaller than those required to 
decouple 3D superconductivity^ 
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C. Coherence length from the VL form factor 

In the case of perfect crystalline VL order, the ob- 
served intensity of the first order diffraction peak 1(H) oc 
X^VIGI where T is the form factor of a single vortex 
and has units of field. The sum is over all the q- vectors 
contributing to the intensity near wave- vector |G|, and 
we have assumed that the rocking curve width remains 
constant with field in obtaining 1(H). A variational so- 
lution to the Ginzburg-Landau model, namely the Clem 
model for the VL form factor j 5 - 2 - yields T cx GK\(G£), 
where K~x denotes the modified Bessel function of first 
order, G = |G| = 2ir^H/a&o, and the vortex core size 
£ is the only fit parameter.^ 3 - We point out that the ap- 
plication of the Clem model to our data yields an up- 
per bound for £; disorder effects are, for example, not 
included.-^ 4 - With increasing vortex lattice disorder, the 
degree by which the Clem model overestimates the coher- 
ence length is larger. By comparing the doping depen- 
dence of the extracted £ with the G inzburg Landau co- 
herence length £gl = \J ^o/2ttH C 2 estimated indirectly 
from specific heat^ and high-field magneto-resistance 
experiments^ a reasonable agreement is found on the 
overdoped side, see Fig. 5. This suggests that our SANS 
data indeed provide a measure of £, even though disorder 
is undoubtedly present. Identifying £ with the Pippard 
coherence length £ p ~ Hvf (k)/A(k) suggests that the rel- 
atively short coherence length £ around optimally doped 
LSCO is not only due to the large pairing gap A; the 
small Fermi velocity is also playing a significant role. 



D. Disorder effects and structure factor 

On the underdoped side, we find a strong discrep- 
ancy between the coherence length estimated from spe- 
cific heat and the SANS data fitted with the Clem model 
for the form factor, see Fig. 5. For LSCO x = 0.145, we 
find £ ~ 80 A corresponding to an unrcalistically small 
upper critical field HqH c2 = $ /27r£ 2 ~5T. A larger co- 
herence length may result from the weakening of super- 
conductivity due to competition with, for example, mag- 
netism. However, this does not explain the discrepancy 
between our SANS data and the specific heat data.— A 
more plausible explanation is that the VL disorder po- 
tential increases with undcrdoping. It is possible that VL 
disorder proliferates as the system is tuned towards the 
state where magnetism and superconductivity coexist. 

The VL displacements throughout the doping range 
are well-described by elastic theory, namely the Bragg 
glass (BrG) paradigm ! 16 ' 17 In the presence of disorder, 
the positional order of an elastic VL decays exponen- 
tially with a characteristic length scale Ra- (With in- 
creasing disorder Ra — > 0.) If such an exponential de- 
cay were to persist at all length scales R, a total de- 
struction of long-range order would result i 5 - 5 - This pro- 
posed destruction, even under weak disorder, presented 
a long-standing puzzle with respect to experimental ob- 



servations where Bragg peaks may readily be observed. 
Theoretically, the puzzle was resolved with the advent of 
the BrG paradigm, in which an asymptotic regime for 
R > Ra enters, where the positional order decays only 
weakly, leading to algebraically diverging Bragg peaks. 
At the crossover scale Ra vortex displacements are com- 
parable to the lattice spacing clq = yj cr<5>o/H. In earlier 
muon spin rotation (/xSR) work on a LSCO x = 0.105 
sample, an order-disorder transition was observed and 
associated with a transition out of the quasi-long-range- 
ordered BrG phased VL correlations in the BrG phase 
have been explored more directly in a study of low-purity 
niobiumi^i We point out that the BrG paradigm can ex- 
plain both the H~ 5 and H~ 2 field-dependences of in- 
tensity (c.f. Fig. 4(b)), as well as the crossover between 
them. Dependent on whether the instrumental resolu- 
tion s « 60ao is larger or smaller than Ra, a different 
field dependence of the SANS intensity is predicted^ 6 - 
For s < Ra, the contribution to the structure factor is 
identical to that of a crystalline VL, hence / ~ \j\fH~. 
In the other limit s > Ra, an additional factor H^^ con- 
tributes to the intensity. Elastic theoryi 6 - yields \i = 3/2 
and hence / m H~ 2 — as indeed we observed previ- 
ously in LSCO x = 0.105^ see Fig. 4b. The intensity 
for x = 0.145 is also consistent with a H~ 5 to H~ 2 
crossover, see Fig. 3. At the crossover field (~ 0.85 T), 
Ra ~ s sa 60ao ~ 3 /jm. By contrast, in LSCO x = 0.105 
the crossover field <C 0.05 T, implying that the disorder 
potential increases dramatically with decreasing doping. 
This disordering seems to occur as the static magnetic 
correlation length ( approaches the VL spacing (£ — > a Q ), 
suggesting an electronic origin to the VL disorder effectu- 
ated by real space competition between £ and clq, rather 
than orthorhombic twin boundaries or impurities as ob- 
served in other superconductors .21 We noticed that the 
BrG model also provides an excellent description for the 
(previously unexplained) field dependence of intensity in 
Nd2_ :c Ce x Cu04 (NCCO) with x = 0.15 (Fig. 4(b)) .21 
More experimental studies are required on NCCO to es- 
tablish if the VL disorder therein has origins similar to 
that in LSCO. 



E. Real space picture 

We now consider how magnetic and superconducting 
order parameters might coexist in real space. Around 
x = 0.12 doping in LSCO, fiSR measurements revealed 
magnetic and superconducting fractions that together ex- 
ceed 100%. It was therefore concluded that the magnetic 
and superconducting order parameters are not phase sep- 
arated but rather intertwined on a nanometer scale i 57 ' 58 
This real space picture is not easily reconciled with neu- 
tron diffraction studies that report a magnetic correla- 
tion length of several hundreds of Angstroms . 25 ' 31 ' 33 ' 34 
One possibility is that the charge of the muon induces 
magnetism in which case the /xSR technique overesti- 
mates the magnetic volume fraction^ Here we showed 
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that the vortex lattice becomes more disordered as the 
vortex interspacing approaches the magnetic correlation 
length. This suggests that the magnetic and supercon- 
ducting order parameters are coupled. How magnetism, 
superconductivity, and vortices are arranged in real space 
when the spin correlation length is larger than the vortex 
spacing (£ > ag) is an interesting question that is difficult 
to address with the SANS technique since no observable 
SANS signal is found in that region of the phase diagram, 
see Fig. 1. Direct imaging techniques 1 ^— are undoubt- 
edly more informative in this regime although it may be 
experimentally challenging to probe magnetism and vor- 
tices simultaneously. 

IV. CONCLUSIONS 

In summary, our studies of the vortex lattice in LSCO 
allow us to draw two main conclusions. First, near op- 
timal doping, and far from the SDW instability, the VL 
structure/orientation and the small superconducting co- 
herence length £ (and hence large upper critical field 
H C 2) may both be rationalized as arising from a vanish- 
ing Fermi velocity due to the van-Hove singularity near 



the zone boundary. Second, we find that the fermio- 
logical picture breaks down as the SDW instability is 
approached by tuning cither the doping or the applied 
magnetic field. There, the vortex lattice structure factor 
needs to be accounted for in the field-dependence of the 
observable SANS intensity. The Bragg glass paradigm, 
describing vortex lattices in the presence of weak dis- 
order, accounts for the SANS intensity behavior across 
the entire phase diagram (where a SANS signal is dis- 
cernible), from the underdoped x =0.105, to the opti- 
mally doped x ~0.16 and overdoped x =0.23 regimes. 
In particular, it is able to explain the cross-over as the 
SDW instability is approached. Evidently, the SDW or- 
der acts as a novel electronic provenience of disorder on 
the vortex lattice in LSCO. 
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dermayer, B.M. Andersen, and N.B. Christensen and 
support from the Swiss NSF (through NCCR, MaNEP, 
and grant Nr 200020-105151, PBEZP2-122855), Dan- 
Scatt, UK EPSRC, and by the Ministry of Education 
and Science of Japan. 
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